FIXED-ENERGY MULTI-PARTICLE MSA 
IMPLIES DYNAMICAL LOCALIZATION 



VICTOR CHULAEVSKY 



Abstract. This work is a continuation of 15 where we described two elemen- 
tary derivations of the variable-energy MSA bounds from their fixed-energy 
counterparts, in the framework of single-particle disordered quantum parti- 
cle systems on graphs with polynomially bounded growth of balls. Here the 
approach of | 15| is extended to multi-particle Anderson Hamiltonians with in- 
teraction; it plays a role similar to that of the Simon— Wolf criterion for single- 
particle Hamiltonians. A simplified, fixed-energy multi-particle MSA scheme 
was developed in our earlier work [7J, based on a multi-particle adaptation 
of techniques from Spencer's paper 1361 . Combined with a simplified variant 
of the Germinet-Klein argument |29l described in 1141 . the outcome of the 
fixed-energy analysis results in an elementary proof of multi-particle dynam- 
ical localization with the decay of eigenfunction correlators faster than any 
power-law. 



1. Introduction. 

This manuscript is an extended version of my talk given recently at the work- 
shop "Mathematics of quantum disordered systems" organized by the Institut de 
Mathematiques de Jussieu at Chevaleret (Paris) and by the Institut Galilee at Vil- 
letaneuse. A reader familiar with |15] may notice that substantial portions of the 
text are borrowed from that preprint, often by the "copy-n-paste" method. A good 
time-saving tool, it is also a well-known source of a.s. creation of notational incon- 
sistencies, for which I sincerely apologize in advance. In my talk, this text has been 
promised to be uploaded promptly to arXiv; a more detailed version will be made 
available later. 

I would like to stress here a point which may have been insufficiently emphasized 
in the talk: although several approaches to the derivation of the spectral and dy- 
namical localization from the FEMSA (fixed-energy multi-scale analysis) estimates 
for single-particle systems are known by nowQ (cf., e.g., [5], [30]), this text focuses 
only on two methods which seem to provide the shortest path to the multi-particle 
localization from the FEMPMSA results, which are much simpler to obtain than 
their variable-energy analogs, the VEMPMSA (variable-energy multi-particle MSA) 
bounds. The choice is certainly biased by the author's personal preferences, and 
there is no doubt that alternative, more efficient and more general methods will 
appear in near future. 



H thank Peter Hislop for a fruitful exchange on this subject during the above mentioned 
workshop " Mathematics of quantum disordered systems" . 
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Several ideas employed in the derivation "FEMPMSA =► VEMPMSA" have been 
used by other researchers in the context of single-particle systems. The "disorder- 
energy" measurable space x K appears already in the work by Martinelli-Scoppola 
[35) . In a more general context, it becomes the main scene of action in the Simon- 
Wolf paper [37] (see [37J for a more extensive bibliography and a discussion of a 
series of works which initiated the research by Simon and Wolf). The fact that the 
matrix elements of resolvents are rational functions with a " moderate" number of 
poles is explicitly used, e.g., by Bourgain-Kenig [5 a , in a "hard" situation (Bernoulli 
Anderson Hamiltonians). Germinet-Klein |30j use spectral reductions for random 
Hamiltonians with virtually no assumption on the marginal probability distribution 
of the random potential (the scatterers amplitude of an alloy potential may follow 
any probability law not concentrated on a single point). While the singular nature 
of the random potential forces one to make use of technically involved analytic and 
probabilistic tools, making a strong assumption on the regularity of the marginal 
distributions (viz., sufficient regularity of the probability density) gives rise to a 
significant simplification of spectral reductions (FEMSA => VEMSA). 

This list can be continued ... 

Needless to say that the Fractional Moment Analysis (FMM) always starts as a 
fixed-energy analysis. 

In the framework of the multi-particle Anderson Hamiltonians, the above men- 
tioned strong assumption on the regularity of the marginal density remains so far 
the only means to achieve more optimal, physically reasonable bounds on long-range 
charge transfer processes (tunneling) in an interacting quantum system of N > 3 
particles; cf. a brief discussion in subsection 12. 31 Resonances occurring in such sys- 
tems can be qualified as " structural" ; they do not appear in 1-particle systems and 
can be treated in a relatively simple way in 2-particle systems. The solution to this 
problem (explicitly analyzed by Aizenman and Warzel [3]), proposed in our works 
[llj . |12j . requires the above mentioned regularity assumption. Therefore, limit- 
ing the spectral reduction (FEMSA => VEMSA) to random Hamiltonians obeying 
this assumption does not seem to be an overly big concession - at least, until a 
new, more efficient solution is found to the problem of long-range tunneling due to 
structural resonances in systems with N > 2 particles. 

On the other hand, proving merely the Anderson localization phenomenon in an 
interacting quantum system with an arbitrary (but fixed) finite number of particles 
is a simpler task than proving efficient decay bounds for eigenfunction correlators. 
This can be done under a much weaker hypothesis of Holder (or even log-Holder) 
continuity of the cumulative marginal probability distribution function (PDF) of the 
external random potential featuring the IID, IAD ^Independence At Distance) or 
an (appropriate) strong mixing property. The main geometrical tool here is the 
notion of " separability" of pairs of finite volumes developed in our joint works with 
Yuri Suhov [5J[9] (for N = 2 particles) and in [10 (for N > 2 particles). 

2. Basic notations, facts and assumptions 

Throughout this paper, we work with discrete Schrodinger operators (DSO) act- 
ing in Hilbert spaces of square-summable complex functions on connected countable 
graphs. Indeed, the techniques and results of the MSA, initially developed for op- 
erators on periodic lattices, are naturally extended to more general graphs with 
polynomially bounded growth of balls (such graphs as Bethe lattices remain so far 
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out of the MSA's reach). Another motivation for presenting the new approach on a 
graph comes from the fact that the natural language for the description of a system 
of N > 1 interacting indistinguishable quantum particles (bosons or fermions) is 
that of a symmetric power of the configuration space Z of the respective single- 
particle system; already in the case where the configuration space is Z = Z d , d > 1, 
its N-th symmetric power is no longer a periodic lattice. 

2.1. Graphs, configurations and graph Laplacians. Consider a finite or count- 
able connected graph (G,£), with the set of vertices Q and the set of edges £; for 
brevity, we will often call Q the graph, omitting the reference to £. We denote 
by dg(-, •) (sometimes simply by d(- , •)) the canonical distance on the graph Q: 
dg (x, y) is the length of the shortest path x ~> y over the edges. We will assume 
that the growth of balls Bl(x) := {y : dg(x,y) < L} is polynomially bounded: 

sup|B L (x)| < C d L d , L>1. (2.1) 

In particular, the coordination number ng(x) := {y : dg(x, y) = 1} of any vertex x 
is bounded by d (even by d — 1). 

Given a connected graph (Z,£z) serving as the configuration space of quantum 
particles, the configuration space of a system of N > 1 distinguishable particles is 
the cartesian product Z N ; it is usually endowed with the following graph structure: 
a pair (x, y) is an edge iff, for some j e [1, N], (xj , yj ) G £z, while for all i ^ j Q 
one has Xi = yi. In other words, y is obtained from x by moving exactly one 
particle to one of its nearest neighbors in Z. To make explicit this choice, we use 
the boldface notations Z N for the vertex set of the A-particle configuration graph 
and £( N ) for the described edge set. In general, boldface notations will be reserved 
for "multi-particle" objects. 

Introduce the following mapping from Z n , n > 1, to the collection of finite 
subsets of Z: 

II : x = (xi, ...,a; n ) !-)■ {xi, . . .,x n }. (2.2) 

We will call LIx the support of the configuration x G Z n . (Note that in the frame- 
work of indistinguishable particles, only LIx would be physically observable.) Simi- 
larly, the support of a "polydisk" B^(x) = x™ =1 Bl(xj ) is the set 

IlB L (x) = U? =1 B L (i j )c2. (2.3) 

Further, given a non-empty index subset J C [1,N] PlZ, define a partial projection 
(or partial support) 

Uj : (x!,...,x N ) i-> {xj,j E J} C Z. 

For J = set, formally, n x = 0. 

Apart from the graph distance d(- , •) = d z n(- , •) on Z N , it will be convenient 
to use the max-distance p and symmetrized max-distance p s defined as follows: 

P(x,y) = max d z (xj,yj), p s (x,y) = min p(x,7r(y)), 

where the elements it e & n of the symmetric group 6 n act on vertices x e Z N 
by permutations of the coordinates. In terms of p, a polydisk B^(x) is a ball of 
radius L centered at x. (If Z = Z d with max-distance, then polydisks are cubes.) 
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The canonical (negative) graph Laplacian (— Ag) on a finite or countable graph 
(G, £) is given by 

(-A g f)(x) = (/(*) - /(»)) = "oW/W - E /(f) ( 2 - 4 ) 

where we use a popular notation (x, y) for a pair of nearest neighbors x,y £ Q, i.e., 
dg(x,y) = 1, and ng(x) is the coordination number of the point x. For brevity, 
we will sometimes use slightly abusive notations like (x, y) G A, A C (? instead of 
(x,y) G (A x A) n Eg. 

From this point on, unless otherwise specified, we will use the notation Q only for 
finite connected graphs, either in the single-particle context or in a situation where 
the nature of a graph (single- or mult i- particle) is irrelevant, while Z will stand 
for a countable connected graph with polynomial growth of balls. Finite connected 
subgraphs of Z N will be sometimes denoted by Q. 

In operator form, we can write, for an arbitrary (connected) graph Q, 

^G ^G ^ ^ ^-x,yi ^x,y I ^-x){^-y \i 

{x,y) 

where ng is the operator of multiplication by the function x i— > ng(x). Given 
a subgraph A C 5, define its internal, external and the so-called edge boundary 
(relative to Q) as follows: 

d g A = {yeA: dg(x, Q\A) = 1}, d+A = d g g\A, 
dgA = {(x,y)ed g Axd+A: dg(x,y) = 1}. 

Working with a given graph Q(c Z), we always mean by a ball B^(it) C Q the 
set {y G : dg(u,y) < R}, i.e., the ball relative to the metric space (Q, dg). 

The Laplacian (hence, a DSO) in a subgraph A C Q can be defined in various 
ways. The two most popular choices are: 

• The canonical (negative) Laplacian in A, (— A\f) — (— Aa/), defined as in 
(|2.4D with Q replaced by A. It this context, it is usually considered as an 
analog of the Neumann Laplacian, and reads as follows: 

(-A^/)(x)=n A (x)- E ^ ( 2 - 5 ) 

{x,y)<£A 

• The Dirichlet Laplacian (-A^ g ) = 1 A (-A^ g ) 1 A \ £ 2 (A). Here we use a 
natural injection £ 2 (A) l 2 (G)- The Dirichlet counterpart of (|2.6I) is 

(-A?/)(x) = ng{x) - £ Hv), (2-6) 

(x,y)EA 

with ng(x) > n A (x), so (— A^) > (— A^) in the sense of quadratic forms. 

We will use the Dirichlet Laplacians and DSO H®. Given a decomposition 
g = A U A c , A c := g \ A c , we can write 

~ Ag = ng - E ^x,y — E ^x,y — E (^x,y+^y,x) 

{x,y)eA {x,y)eA c (x,y)£dA 

= ((_A°)©(-A°))-r A , s 
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with r A , e = T,( XlV )edA ( r *^ + r v^)- Respectively for the DSO Hg = -Ag + V, 
where V : Q — > K is usually referred to as the potential, one has 

Hg = H' gA - r A , g , Hg,A ■= ("Ag + V) © + V). 

We omit the superscript " N" , since the nature of the boundary conditions in Q is 
not related to the choice of Dirichlet or Neumann decoupling induced by Q = AuA c . 

The spectrum of a (finite-dimensional) operator Hg, i.e., the set of its eigenval- 
ues (EVs) counting multiplicities, will be denoted by T,(Hg). The resolvent of a 
Hamiltonian Ha(w), A C Z N , will be denoted by Ga(E;w) and its matrix ele- 
ments in the delta-basis (Green functions) by Ga(x, y; E; u). The subscript A will 
be omitted when A = Z N . 

In a number of formulae and statements, we will use parameters j3, r, g £ (0, 1), 
and a £ (1,2). Unless otherwise specified, we assume that f3 = 1/2, r = 1/8, 
g = (a — 1)/2 = 1/6 and a = 3/2. Note that the exponent figuring in Definition 
l2~T1then equals 7/8. 

The function 7 : (to, L) <-¥ to(1 + L~ T ) introduced below is a convenient replace- 
ment for the decay exponent ("mass") used in the MSA and dependent upon the 
scale. We merely make this dependence explicit. Clearly, 7(771, L) > m for any 
L>0. 

Definition 2.1. Given numbers E e K, to > and L e N*, a ball B/_,(u) is called 

• _E-resonant (-E-R, in short), if dist(E(HB i ( u )) 5 E) <e~ L , and _E-nonresonant 
(E-NR), otherwise; 

• (E, TO,)-nonsingular ((E, m)-NS), if for all x, y 6 Bj,(u) with p(x, y) > L~^~ 

|aBi(u)| • |G Bjl („)(x,y;£OI < e^™^"^, (2.7) 

where 

7 (to,L) :^TO(l + i- T ), (2.8) 
and (E, m)-nonsingular ((E, m)-NS), otherwise. 

2.2. Geometric resolvent inequality. The second resolvent identity implies the 
so-called Geometric resolvent equation for the resolvents Gg(E) = (Hg — E)^ 1 , 
G A ,(E) = (H A < - E)~\ Gg A (E) = (H' A - E)~ l : 

Gg{E) = G'g(E) + G'g{E)T Kg Gg(E). (2.9) 

For x,u e A and y € A c , one has Gg(x,u;E) — G^(x,u;E) and G g (x,y;E) = 0. 
This results in the Geometric resolvent equation for the Green functions 

Gg(x,y;E)= £ G%(x, u; E) Gg(u', y; E) (2.10) 

(u,u')EdgA 

and the Geometric resolvent inequality (GRI) 

\Gg(x,y;E)\< £ \G%(x,u;E)\ \Gg(u',y;E)\. (2.11) 

In the rest of the paper, the superscript "D" will be omitted, since we do not use 
the Neumann boundary conditions. 
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2.3. Assumptions on the random potential. The efficiency of the quantitative 
bounds on eigenfunctions and eigenfunction correlators for N > 3 interacting par- 
ticles available at the moment and the complexity of the proofs depend upon the 
assumptions on the random potential. 

For clarity of presentation, we always assume that the random potential field 
V : Z xQ^Rona graph Z is IID. 

The first, more general condition (leading to weaker results) is as follows: 

(Wl): The marginal probability distribution function (PDF) Fy is uniformly 
Holder continuous: there are constants Ch £ (0, +oo), S S (0,1] such that 

sup(F v (t + s) - F v (t)) < C H s S . (2.12) 

tea 

To obtain more optimal decay bounds on eigenfunction correlators, with the help 
of a method developed in [T2]-[T1], one needs an additional assumption which we 
will describe now. Introduce the following notations. Given a finite subset A C Z, 
let £a(w) be the sample mean of the random field V over the A, 

£ A (u) = \A\- 1 J2v(x,u } ) 

xeA 

and define the "fluctuations" of V relative to the sample mean, n x (uj) = V(x,ui) — 
£a(w), igA. Denote by the sigma-algebra generated by {rj x : x G A}, and by 
-Ff A (- \3t], a) the conditional distribution function of £a given 3^, A- Assuming that 
AcZ and diam(A) < R, introduce, for s > 0, an g^A-measurable random variable 

i>£ a (s;uj) := ess sup sup \F^ A (t + s |&,,a) - ^ A (t |&,,a)|- (2-13) 

The stronger assumption on the random field V is as follows: 

(W2): There exist C ,C" , A' , A" , B' , B" e (0, +oo) such that for any finite 
subset A C Z with diam(A) < R, 

Vse(0,l] ¥{u $a (s;uj) > CR A s B ) < C R A ' s B ' . (2.14) 

In the case where the random field V is IID, a more natural condition should 
refer to the cardinality A| of the subset A, but we will use the above property only 
in a situation where A is a subset of some ball, of an explicitly bounded diameter. 

For further use, define the following function on [0, 1]: 

h L = h[ N ' x ' y) : s i — y |Bf } (x)| • \B[ N) (y)\CL A s B + C'L A 's B '. (2.15) 
See the discussion of the validity of (W2) in subsection 18.21 below. 

2.4. Assumptions on the interaction potential. For brevity, we consider only 
finite-range interactions U generated by a pair interaction potential U^(-): 

(Ul): There is an integer r$ > and a function : N — > K with finite 

support, suppC/^ 2 -* C [0, 7*o], such that 

U(x)= UW(d( Xi ,Xj)). (2.16) 

l<i<j<N 

The parameter rp will be called the range of the interaction U. 
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In the case where the single-particle configuration space is a group (e.g., Z d or 
R d ), the above form of U corresponds to translation invariant interactions; in the 
case of a Euclidean space, such an interaction is also isotropic. In all cases, U of the 
form (|2.16|) is permutation symmetric. We stress that neither of these properties 
is crucial for the multi-particle MSA scheme. Moreover, this scheme can be easily 
adapted to interactions with a hard core. 

The assumption of finite range simplifies the induction on particles. In |13| . we 
described an extension of the variable-energy MPMSA to interactions decaying at 
an exponential or subexponential rate, 

|C/ (2) WI < cie" C2rC3 

with c±,C2 G (0, +00) and C3 > sufficiently close to 1. The proof of the fixed- 
energy MPMSA bounds can follow essentially the same path. In fact, the most 
significant modification is required in the proof of Lemma 16.41 see the details in 
[T5] . In the framework of the fixed-energy MPMSA (which is simpler) , we plan to 
address infinite-range interactions in a forthcoming paper. 

3. Main results 
Let 38 \ be the set of all bounded Borel functions / with ||/||oo < 1. 

Theorem 3.1. Fix any integer N > 2. Under the assumptions (Wl) and (Ul), 

for any N G [1, N] there exists go(N) < +00 such that for all \g\ > go{N) and some 
m = m(g, N) > Const(N) In \g\ > 0: 

(A) with probability one, the random operator H(cj) = Ho + g'V(uj) + U has pure 
point spectrum and all eigen] 'unctions *&j(u>) rapidly decaying at infinity: for each 

j and some Xj, for all x and some C(x.j,ui) G (0, +00), 

l*^)! < C(Si j ,u;)e- al « 1+c r^; (3.1) 

(B) for all points x,y and some a,c, C(x) G (0, +00), 



E 



sup |{l y |/(H(w))| 1, 



< C7(x)e- alnl+Cp ^ y ). (3.2) 



With f — fx : t h-> e ltx , A G R, Eqn (|3.2[) reads as the strong dynamical 
localization for the random Hamiltonians H( n \lj). 

Technically, we prove upper bounds in p.l[) - (|3.2l) of the form LI K ( 1+e ) at dis- 
tances ~ Lk ~ (Lo) a , which can be easily translated into the above decay rate. 

Assertion (A) can be made stronger, without replacing the hypothesis (Wl) by 
a stronger condition (W2). Indeed, one can prove an exponential decay of all eigen- 
functions, using the variable-energy MPMSA developed in [TU]. In the present 
paper, we focus on a simpler approach, based upon the fixed-energy MPMSA. 

Note also that the decay bounds I3.1H3.2I have an unusual form as compared 
to their more traditional counterparts from the single-particle MSA: constants 
C(xj, w), C(x) depend upon the multi-particle configuration x (actually, upon the 
diameter of its support, diam({xi, . . . , Xjv}))- This feature is present in the first 
papers on multi-particle localization |10j . [3]; we discuss it below in more detail. 

As was already said, in order to obtain more optimal localization bounds, viz.: 
exponential spectral localization (as in [TU]) and uniform estimates on the decay of 
eigenfunctions and their correlators, we make use of a stronger assumption (W2). 
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Theorem 3.2. Fix any integer N > 2. Under the assumptions (W2) and (Ul), 

for any N G [1, N] there exists go(N) < +oo such that for all \g\ > go(N) and some 
m = m(g, N) > Const(N) In \g\ > 0: 

(A) with probability one, the random operator H(cj) = Ho + gV(w) + U has pure 
point spectrum and all eigenf unctions ^fj(uj) decay exponentially fast at infinity: for 
each ^fj, some Xj(aj) G Z N , Cj(u>) € (0, +oo) and all x G Z N , 



|*i(x,w)| < Cj(w)e 
(B) for all points x, y and some a, c, C G (0, +oo) 



-mp s (x,Xj). 



IE 



sup |(l y |/(H(w))| 1, 



< C( 



(3.3) 



(3.4) 



It is readily seen that, given a point u G Z N (e.g., u = Xj), 

P(X, U) > p(x, {7r(u), 7T G ©tv}) — > +oo, 

p(x,u)— f+co 

so that the assertions (A) and (B) of Theorem 13.21 imply their counterparts from 
Theorem 13.11 (under a more restrictive hypothesis (W2)). An important point is 
that Theorem 13.21 provides a more precise and explicit information on the decay 
properties of the EFs and their correlators. Moreover, once formulated in a physi- 
cally relevant language of indistinguishable particles (bosons or fermions), Theorem 
13.21 operates with a natural distance p s between particle configurations. Possible 
lack of decay in terms of the non-symmetrized distance p is due to the fact that 
any zb-symmetric eigenfunction \& (x) takes the values of the same magnitude along 
an orbit <^(x) := {n(x), n G &n}- 

Working with restrictions of the Hamiltonian H^^u;) to the subspaces of sym- 
metric or antisymmetric functions on Z N would be simpler on the A^-th symmetric 
power of the 1-particle (physical) configuration space Z. In the framework of the 
variable-energy scaling analysis, such an approach has been used in our earlier work 
|13j . Here, for the sake of brevity, we prefer to work with distinguishable particles, 
i.e., with Hamiltonians H^ n '(lj) in the entire Hilbert space £ 2 (Z N ). This also leads, 
formally, to more general results. Nevertheless, it is to be stressed again that only 
the language of indistinguishable particles is physically relevant. 

Finally, note that adaptations of the variable-energy MPMS A to localization near 
the lower edge of the spectrum (without the assumption of strong disorder) and to 
weakly interacting systems ("sufficiently" localized without interaction) have been 
obtained by Ekanga [2"Tll22) . in the case where Z = Z d , under the assumption of 
log-Holder continuity of the PDF Fy of the external random potential. 

Theorems 13 . 1 1 and 13.21 give the final results of the multi-particle localization anal- 
ysis. However, the novelty of the present work resides essentially in Theorems 19.21 
and 19.31 deriving the spectral and dynamical localization from the results of the 
fixed-energy analysis of resolvents. 

4. SUBHARMONICITY ON GRAPHS 

The results of this section apply to arbitrary graphs, including Z N , N > 2. 
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Definition 4.1. Let Q be a finite connected graph, L > £ > two integers and 
q G (0, 1). A function / : Q — >• R + is called (£, g)-subharmonic in a ball B/,(u) C 
if for any ball B^(x) C Bj(u) one has 

/(z)<? max f(y). (4.1) 

We will often use the notation A4(f, A) := max xe A |/(^)|- 

Lemma 4.1. //a function f : Q —> R + defined on a finite connected graph Q is 
(£,q)-subharmonic in a ball C with L > I > 0, then 

f(x)<ql^\M(f,g)<q^M(f,g). (4.2) 

In fact, the factor M(f, Q) in the RHS of (|4.2p can be replaced by M(f, B L+1 (x)). 

Proof. See [H]. □ 

Lemma 4.2. Consider a finite connected graph Q and a ball B/,(u) C Q 7 with 
L > £ > 0. Fix numbers E G R, m > and suppose that all balls Re(x) inside 
Bi(u) are {E, m)-NS. Then V t/ G Q \ Bl(m) i/ie function 

f i x i ^ IGgf^y;^)! 

is (£, q)-subharmonic in Bi(w) imi/i g = e~ 7 ^ m ^'^. 

Proof. The claim follows directly from the Definition 14. II □ 

Lemma 14.11 suffices to assess the Green functions in a ball B^(u) which does 
not contain any singular £-ball, but to analyze the situation where B^(u) does not 
contain any pair of disjoint singular £-balls, one needs the following extension of 
Lemma 14.11 inspired by the proof of Theorem 1 in [55] : approaching a single " bad" 
ball separately from the points x and y. 

Lemma 4.3. Let Q be a finite connected graph, and f : Q x Q — > R + , / : (x,y) i— > 
f(x,y), be a function which is separately (£,q)-subharmonic in x G B r '(u') C Q and 
iny G B r » (it") C 0, with r',r" >£>0 and d(u', u") > r' + r" + 2. TTiera 

/(u',u") < ql^\ + [^\M(f,g x 0) < /^""^(/.gx 0). (4.3) 
Proof. See [HI. □ 

5. Eigenvalue concentration bounds 

5.1. One- volume EVC bounds. The first eigenvalue concentration (EVC) bound 
required for our scheme is a direct analog of the Wegner estimate 40] for multi- 
particle operators. It can be proved under the assumption (Wl). 

Theorem 5.1. Let V : Z x f2 — > R be a random field satisfying the hypothesis 
(Wl). Then for any (3' G (0,/?) and large Lq, the following bound holds true: 

VxeZ" sup P{ B L (x) is E-R }<e- Lli ' (5.1) 



With (3 = 1/2, it suffices to set /?' = 1/4. 

Naturally, Theorem 15.11 is based on a probabilistic bound on the events of the 
form {lo : dist(_E, S(Hb l ( x ))) < s }, with s = nT^ . An optimal Wegner-type 



10 



V. CHULAEVSKY 



bound, for an IID random potential V with Lipshitz-continuous marginl PDF Fy, 
would have the form 

P { lo : dist(E, E(H Bt ( x ))) < s } < Const|B L (x) s. (5.2) 

Indeed, such a (multi-particle) bound was proved by the author and Yuri Suhov 
for IID random potentials with PDF analytic in a strip {E G C : | Im £7 1 < e}, 
for some e > (hence, having an analytic probability density). This covers the 
so-called a-stable laws indexed by a real parameter a e [1,2] (not to be confused 
with the scaling exponent a appearing in the recursion Lk+i — [L%]), where a = 2 
corresponds to the Gaussian law and a = 1 to the Cauchy law. The proof is 
based on Molchanov's rigorous path integral formula for the unitary propagators 
e -itH (j QT finite-difference operators H) and follows essentially the same path as 
in Carmona's proof of analyticity of the Density of States (DoS) for tight-binding 
Anderson Hamiltonians; cf. |18j . The same method applies to the two- volume 
bounds given in Theorem 15.31 below, again for analytic PDF Fy. 

Kirschj [33] proved an optimal multi-particle Wegner estimate of the form (|5.2[) 
for uniformly Lipshitz continuous marginal probability distributions. The proof 
of Theorem 15.11 in 8 , for possibly singular (but uniformly continuous) marginal 
probability distributions is based on Stollmann's lemma from [35] which gives rise to 
a non-optimal volume dependence; the latter, however, is absorbed in the fractional 
exponent and is more than sufficient for applications to the multi-particle MSA. 
Stollmann's lemma itself is very general and not sensitive to the regularity properties 
of the marginal distribution, but is effective (suitable for MSA applications) only for 
marginal PDFs featuring at least log-Holder continuity, with suitable parameters. 

A one- volume EVC multi-particle bound for continuous systems (Anderson Hamil- 
tonians with an alloy- type random potential), with optimal volume dependence, has 
been reported by Hislop and Klopp |31] . A variant of such bound (again, for alloy- 
type potentials in R d ) based on Stollmann's lemma, with a non-optimal volume 
dependence, was proven earlier in [B] and used in |16) . in the proof of multi-particle 
dynamical localization for Anderson- type Hamiltonians in R d . 

5.2. "Separability" in the multi-particle configuration space. All above 
mentioned one-volume EVC bounds suffice for the proof of the a.s. exponential 
decay (hence, square summability) of multi-particle Green functions. This point 
was not clear earlier, since Simon- Wolf method from [37] does not apply to multi- 
particle Hamiltonians. However, both the variable-energy MPMSA from JTU] and 
the derivation of the VEMPMSA bounds from their fixed-energy counterparts (pre- 
sented here) require two-volume EVC bounds which are more sophisticated for 
multi-particle systems than in the conventional, single-particle localization theory. 
It turns out that the positions of the balls play an important role. For this reason, 
following [10:, we introduce 

Definition 5.1. A ball B^(x) is separable from a ball B^(y) if there exists a 
decomposition B^°(x) = B^' } (x') x B^"" ) (x") such that (cf. Eqn (|23])1 

nB^' } (x') n nB^ ( y ) = 0. (5.3) 

A pair (B^- 1 (x) , B^ (y)) is separable if one of the balls is separable from the other. 



2 See also the preprint arXiv:math-ph/0704.2664 (2007). 



FIXED-ENERGY MULTI-PARTICLE MSA IMPLIES DYNAMICAL LOCALIZATION 11 



Pictorially, (|5.3p says that there is a sub-sample of the potential which determines 
the operator H Bi ( y ) completely, while H Bl ( x ) (more precisely, its component rela- 
tive to a sub-configuration x') has a non-degenerate (and controllable) probability 
distribution, conditional on Ha^fy)- 

Lemma 5.2. For any x G Z N there is i?(x) such that any ball B^^y) with 
p( x i y) > * s separable from B^(x). 

Proof. Let r — diam(IIx) and set R = i?(x) = r + 3NL. Consider a ball B^^y) 
with p(x, y) > R. By definition of the max-distance, there is at least one particle 
position yi G ITy such that d(yi ,Xi ) > R. Consider the maximal connected 
component A y := U^j^hilli) of the union U^E^l^) containing yi a ; its diameter 
is bounded by 2NL, and by triangle inequality, 

dist (A y , ILB L (x)) > R - (diam ILx + 2L) - diam A y > 0. 

Taking as the subconfiguration y' the union Ui<zj{yi}, we conclude that B^(y) 
is indeed separable from B^ (x). □ 

Theorem 5.3 (Cf. Lemma 2 in [E]). Under the hypothesis (Wl), the following 
bound holds true for any pair of separable balls B^(x), B^(y) and all s G [0, 1]: 



Lemma 2 in [TU] is formally stated for operators on a lattice (l* d ) N , but its proof, 
based on a very general Stollmann's lemma on monotone functions }38U39j . applies, 
with minor notational modifications, also to Hamiltonians on graphs, provided that 
the extrenal random potential field V : Z x fi — > K is IID with Holder-continuous 
marginal distribution. 

The dependence of the upper bound on the EF correlators in Theorem 13.11 upon 
the positions x 3 , x is explained by the nature of the two- volume EVC bound in 
Theorem 15.31 under the hypothesis (Wl), "tunneling" between quantum states x, 
y is ruled out with high probability only if the balls B^(x), B/,(y), with a suitable 
L = 0(p(x, y)), are separable. 

5.3. "Weak separability" at large distances. 

Definition 5.2. A ball B^x) C Z N is weakly separable from B^ (y) if there 
exists a ball A C Z in the 1-particle configuration space, of diameter R < 2NL, 
and subsets J\,Ji C [1,N] (1Z such that |j7i| > | J?2 1 (possibly, J 2 — 0) and 



A pair of balls (B^ (x), B^ (y)) is weakly separable if at least one of the balls is 
weakly separable from the other. 

The physical meaning of the weak separability is that in a certain region of the 
one-particle configuration space, there are more particles from configuration x than 
from y. As a result, some local fluctuations of the random potential V(-;uj) have a 
stronger influence on EVs relative to B^(x) than on EVs relative to B^(y). It is 
easy to see that this condition is indeed weaker than the separability. 



P 



{dist (s(H Bi(x) ),E(H Bi(y) )) <s} 



< Const L {2N+ ^ d s s . 



(5.4) 




(5.5) 
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Lemma 5.4 (Cf. Lemma 2.3 in 11 ). If p s (x, y) > 3NL, then balls B^(x) and 
B^ ^ (y) are weakly separable. 

The lower bound by 3NL is not sharp, but slightly less cumbersome than the one 
established in [IT]; the actual value of the constant (in front of NL) is irrelevant. The 
proof is given in [11], formally, applies to the case where Z — Z d , but the extension 
to more general graphs is straightforward. The main idea can be summarized as 
follows: two balls of radius L with centers x, y are not weakly separable iff all 
particles from x occupy approximately the same positions as the particles from y 
(up to a permutation: this is why the symmetrized distance p s is used) . Specifically, 
there has to be a permutation n 6 &n such that Vi 6 [1 , iV] , d^, = O(L), 

yielding p s (x,y) = 0(NL). 

We will call a pair of balls B^ V) (x), B^ V) (y) distant iff p s (x, y) > 3NL. 

Theorem 5.5 (Cf. Lemma 3.1 in |llj). Under the hypothesis (W2), for any pair of 
weakly separable (e.g., distant) balls Bj,(x), B^(y) the following bound holds true: 

P{dist (£(H Bi ( x )),£(H BL ( y ))) <s} <h L (2s), 

with /il(s) defined in (|2.15|) and constants C,C',A,A',B,B' defined in (W2): 

h L (s) = |B L (x)| • \B L (y)\CL A s B + C'L A 's B '. (5.6) 

The proof given in [11] in the case of a periodic lattice Z — X d extends with no 
difficulty to more general graphs. 

6. Fixed-energy multi-particle scale induction 

It is to be emphasized that the entire fixed-energy scaling procedure, for any 
number N > 2 particles, does not reveal the difficulty, mentioned in the Introduction 
and related to the "structural" resonances. As a result, the exponential decay 
(hence, square-summability) of Green functions at a fixed energy can be established 
under a weaker hypothesis of Holder (or log- Holder) continuity of the marginal PDF 
Fy of the random potential field V : Z x — >• R. 

6.1. Scaling of Green functions in absence of tunneling. 

Definition 6.1. A ball Bi fe+1 (u) is called i?-tunneling (E-T) if it contains two 
disjoint (_E, m)-S balls of radius Lk, and £?-non-tunneling (_E-NT), otherwise. 

Lemma 6.1. // a ball B Lfc+1 (u) is £-NR and E-NT, then it is (E, m)-NS. 

Proof. See the proof of Lemma 5.1 from [T3]. Its statement is deterministic and not 
specific to the analytic nature of the potential, thus applies to DSO on graphs with 
single- and multi-particle structure of the potential. □ 

6.2. Main inductive bound. Given a number k > ^z^Nd, pick any < 9 < 
T^a — ~~~~~~ and introduce the following double sequence: 

P(N, k) = P(N, k; N, k, 9) = 3"- n k(1 + 9) k , N = 1, . . . ,N; k > 0, (6.1) 

and a statement, or property, depending upon integer parameters N £ [1, N], k > 1, 
and an interval ICR, which we are going to prove by induction: 

S(I,N,k): VE eJ Vn e [1,N] Vxe2™ 



FIXED-ENERGY MULTI-PARTICLE MSA IMPLIES DYNAMICAL LOCALIZATION 13 



p{b£J(x) is (E,m)-S } < Lf*~ n « = L~ P ^ k \ (6.2) 

6.3. Initial scale bounds. First, consider the case of strongly disordered systems. 

Lemma 6.2. For any Lq > 2, m > 1, N > 1 and N e [1,N] there is g* — 
<7*(m,p, n, Lq) < oo such that if \g\ > g* , then the property S(I,N,0) holds true. 

The proof is similar to its well-known 1-particle counterpart (cf., e.g., |19| ) and 
will be omitted. It also follows easily from the Wegner-type estimate (cf. Theorem 
I5.ip combined with the Combes-Thomas estimate. 

For energies close to the lower edge of the spectrum, an initial scale bound 
suitable for the purposes of the MPMS A has been proved by Ekanga • 

In the framework of multi-particle systems with interaction, a third scenario 
appears where the localization can be proved, merely by an adaptation of the initial 
scale estimate: systems localized without interaction and perturbed by a sufficiently 
weak interaction (e.g., of short range). Such scenario has been first analyzed by 
Aizenman and Warzel [3] with the help of the MPFMM. Ekanga [H] adapted the 
method of [TU] to weakly interacting systems. In particular, he proved that N- 
particle systems in Z , with two-parameter Hamiltonian 

H = — A + .gV(x; w) + ftU(x) 

are localized for any nonzero g £ R and all sufficiently small \h\ < h*(\g\). Note, 
however, that the random potential is required to have log-Holder continuous mar- 
ginal PDF Fy (with suitable parameters), so singular (e.g., Bernoulli) potentials are 
not allowed in [22]. Indeed, the proof follows the MSA scheme (which is perturba- 
tive, by its nature), for the specifically one-dimensional methods based on Fursten- 
berg theory do not apply even to 2-particle Hamiltonians. The same difficulty is 
encountered in [3H3] where the fractional moment method is also perturbative. 

Lemma 6.3 (Cf. |211I22| ). For any L > 2, m > 1, N > 1 and N G [1,N] there 
is g* = g*(m,p,n, Lq) < oo such that if \g\ > g* , then the property S(I,N,0) holds 
true. 

Proof. (We only sketch the proof.) Owing to the positivity of interaction U, of the 
external potential V and of the graph Laplacian —A, one has 

Therefore, a lower bound for the spectrum of HW is provided, e.g., by the lowest 
eigenvalue of the single-particle Hamiltonian. For the latter, it is well-known and 
is due to the Lifshitz tails asymptotics. □ 

The case of the random potential bounded from below by an arbitrary quantity 
> — oo can be treated similarly, replacing V(x; u>) by V(x; ui) — Eq and by 
jjW — NEq 1. In the case where Eq = — oo (e.g., Gaussian potentials) the proof is 
even easier and follows virtually the same path as for strongly disordered systems. 

For brevity, below we treat only the case of strong disorder; an adaptation 
to the band-edge localization merely requires the control of the "mass" m > and 
simple modifications of scaling formulae. 
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6.4. Green functions in decoupled systems. 

Definition 6.2. An iV-particle ball Bl(u) is called partially interactive (PI) if 
diam(riu) > UNL, and fully interactive (FI), otherwise. 

If the interaction has finite range tq and L > 8ro, a polydisk B^(u) can be 
represented as follows: 

B L (u) = Bf >(u) = BfV) x B?'\u") (6.3) 

where n' + n" = N and the subconfigurations u = (u',u"), u' 6 X n , u" G X n 
satisfy the condition 

p(nBi"V),nBf V)) >r . 



As a result, the operator H^l. in B^(u) is algebraically decomposable in the 
following way: 

Hff (u) =<;) u0 ®l(™'0 + l(»0 8H ^ u „ ) (6.4) 

since, for any x = (x', x") G B^ (u), the interaction between the sub-configurations 
x' and x" vanishes. Such a decomposition, if it exists, may be non-unique. We will 
assume that one such decomposition (referred to as the canonical one) is associated 
in some way with every A^-particle PI ball. 

Lemma 6.4. Assume that the interaction U satisfies the condition (Ul). Fix 
an energy E G M. Consider a PI N -particle ball with canonical decomposition 
B^P(u) = B^ (u') x B^™ ^(u") and a sample of the potential V such that 

(a) B^°(u) is E-NR. 

(b) VA" G E(H_ <„«>, „,) the ball B^'W) is (E - A",m)-NS; 

B Lfc (u ) k 

(c) VA' G S(H B („/) ) the ball B { [p(u") is (E - A',m)-NS. 
Then the ball B^°(u) is (E,m)-NS. 

Proof. See subsection llO.il □ 

Corollary 6.5. Assume the property S(I,N,k) and let Lq > 4C^. Then for all 
E el and any PI ball of the form B^ } (x) C Z N , 

p{ B W(x) is (E,m)-S } < ^IT**™ (6.5) 

Proof. Denote by Sk the event in the LHS of (|6.5I) . Consider the canonical decom- 
position B^(x) = B^x') x B^'^x"). By virtue of Lemma ED 

P { S k } < P | (x) is E-R | + P | BfJ (x) is E-NR and (E, m)-S } . (6.6) 

The first term in the RHS is bounded with the help of Theorem l5.ll so we focus on 
the second term. Apply Lemma 16.41 since the option (a) is ruled out, it remains to 
assess the probability of events listed in options (b) and (c). Consider the former: 

P{5'}^p| 3A"€E(H B £ Vo ): B^(x') is (E - X",m)-S | 

E P< 3 A" G E(H B („")^ x „p : B^ } (x') is (E - A", m)-S ] 
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where J" is the sigma-algebra generated by the values of the random potential in 
B^™ (x"). By definition of the canonical decomposition of a PI ball, 

nB£' ) (x')nnB^" ) (x")-0 ! (6.7) 

and since the random field V is IID, for any E" 6 R, including E — A", one has 

(a.s.) p{ B^'V) is (E",m)-S | 5" } - p{ B^'V) is (E",m)-S } . (6.8) 
On the other hand, by the assumed property S(I, N — 1, k), for n' < N — 1, 

P{ B^V) is (E",m)-S } < C- 2N L- p{N ^ k) = C~ 2N L~ 3P{N > k) . (6.9) 
Therefore, we obtain, with L > 4C^\ 

P{5'} < |b£'V)I sup p{b£V) is (£",m)-S} 

< C d L Nd L -3P(N,k) < i C -2N L -±tP(N,k)-($P(N,k)-Nd-l) ^ ^ 
Cj N 4 

since p > 6iVd, so \P(N, k) > ^Nd >Nd+l (for N > 2). Similarly, 

3 A' e S(H B( „ 0rx ,J: B^'V) is (E-X',m)s) < W™^^. 



(6.11) 

Collecting ([53]) , (jol)]) . (pHD) and ([6"TT]) . the assertion follows. □ 

6.5. Scale induction. Introduce the following notations: with k > 0, 
P fe = sup P{B Lfc (x) is (E,m)-S}, 

Q fe+ i = 4 sup P { B Lk+1 (x) is E-R } , 

Sfc+i = sup P { BL k (x) contains a PI (E, m)-S ball of radius Lk } ■ 
For future use, note that by Corollary 16.51 with k > QNd and < 9 < 1/6, 

c ^ 1 r~T P W fc ) n N T Nd / 1 r -^-§P(JV,fc)-ATd 

2Cf^ k ' d k+1 ACj N k+1 

J_ r "6{T+fl) K ( 1 + e ) ^ 1 r 8(1+^) V ; (ft 191 

1 r -^«(i+9) fc + 1 -JVd 1 K (i+e) fc +! 

— 4C2Af ^fc+l — 4C27V 

Theorem 6.6. Let a e (1,2), k > 2JM* and e (0, 2=3. - 2J¥M) ; assMme ^ e 
condition (Wl). If there is an integer Lq > 4(7^ smc/i i/iai 5(7,^,0) is fulfilled and 
Qo < C~ 2N Lo K , then S{I,N,k) holds true for all k>0. 

Proof. It suffices to derive S(I, AT, fc + 1) from S(I, AT, fc), so assume the latter. By 
virtue of Lemma [6.11 if a ball Bi fc+1 (u) is (E,m)-S, then it is either E-R or _E-T. 
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By Eqn (|6.12p . the probability to have at least one PI (E, m)-S ball of radius 
inside B Lk+1 (u) obeys S fc+ i < j^w L k+i +9) ■ 

Further, there are < \C^ N L^ff pairs of disjoint L^-balls in Bx, Jt+1 (tt), thus 



Pk.+l < TrCj^^fe+lPfc + 7Q*:+1 + Sfc+1- 



By Theorem EH Q fc+1 < Const L^e'^+K The function 

/ : L i-> In (const L~ K ) - In (const L d c~ L ") = L - const In L 



on [l,+oo) is either non-negative or admits a unique zero. In either case, the 
assumption Q < C^ 2N L^ K implies Q fc+ i < Cj 2N for all k > 0, so 

-<4k+l + bfe+i < -O d L k+1 

Therefore, 



n ^ V 2 « T 2Nd v> 2 , V-2JV r -Kii+wj 1 ^ v d" r~ 



^ r -2JV r -K(l+«) t+1 

- ^fc+1 ' 

provided that ^ + 27Vd > k (i.e., /c > ^f ) and = (f - ^) - 1 > 0. □ 

This marks the end of the fixed-energy multi-particle multi-scale analysis. 

7. From fixed to variable energy: First approach 

Now we establish a fairly general relation between fixed-energy probabilistic es- 
timates on the Green functions and variable-energy bounds for two disjoint finite 
volumes. It does not matter how the probabilistic input is obtained; in particular, 
the results of this section can be combined both with the MSA, performed for each 
fixed energy E in a given interval /Ct, and with the FMM (which always starts 
as a fixed-energy analysis). 

7.1. Derivation of the first variable-energy bound. The results of this sub- 
sectioi{| are based on a straightforward adaptation of the techniques developed by 
Elgart et al. [23]. In fact, the assertion of Theorem 17. II below is a mere encapsula- 
tion of an argument from |23| in a statement involving four parameters which can 
be adapted in various ways to particular models. 

It is convenient to assume that |/| = 1, so the interval / with the Lebesgue 
measure mes(-) is a probability space, and so is the product space (f!xJ,Px mes). 
Given L £ N and points x, y G Z N , set for brevity 

F x , y (E) = |GB L (x,y;£)|, F x (£) = max F x , y (£), (7.1) 

yea-Bi(x) 

and introduce the subsets of I parameterized by a > 0: 

<f x , y (a) = {E e I : F x , y (£) > a}, fijp) = {E e I : F X (E) > a}. (7.2) 
(The L-dependence will be often omitted for brevity.) 



^1 thank Tom Spencer and Sasha Sodin for stimulating discussions of the works I23II24| . and 
Ivan Veselic for a fruitful exchange on this subject. 
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Theorem 7.1. Let L > 0, x 6 .Z , y G <9~Bl(x). Let { } x &e the eigenvalues 
of the operator HB t ( x )(w) and I d R an interval of unit length. Let be given 
numbers a,b,c, Vl > such that 

b < min{|B L (x)|- 1 ac 2 , c} , (7.3) 

and for all E £ / 

P{# x (a)}=P{F x (S) >a}<V L . (7.4) 
There is an event B x (b) wif/i P { £ x (6) } < fc" 1 ^ swc/i i/mi V w S x (fe), i/ie set 

<f x (2a) = <? x (2a;w) = {.E : F X (E) > 2a\ 
is contained in a union of intervals UjL-^Ij, Lj := {E : \E — Xj\ < c}, Xj E I . 
Proof. Consider the following events parameterized by b > 0: 

BJp) = {uj e ft : mes(<? x (a)) > b}. (7.5) 
Apply Chebyshev's inequality and the Fubini theorem combined with (|7.4p : 
P{£ x (6)} < ^ 1 E[mes(^ x (a))] 

= b-' J^dEE[ l {Fx( £)>a} ] < b- X V{L). (7.6) 

Now fix any ui B x (b), so that mes(S' x (a);uj) < b. There is a subset {Xj}^^ of 
the EVs of the operator H B£ / X \ such that the Green function E n- G Bl/ ( x )(x, y; E) 
reads as a rational function (below we remove the vanishing terms, if any) 

N' 

f:E^G BUx) ( X ,y;E)=:^2-^, N' < N := |B L (x)|; (7.7) 

j= i i 

here Kj = Kj(x,y) f and £\ |kj-| < £\ |V>i(x)V>i(y)| < A. Let 
M{2c) = {AeM: min |A 3 - A| > 2c}, 
^(c) = {Ael: min |Aj - A| > c}, c> 0. 

Observe that, with < 6 < c, ^4 := {£ : dist(£, M(2c)) < b} C ^(c), hence, the 
set s^h is a union of open sub-intervals at distance > c from the spectrum, and on 
each sub-interval one has \ f'(E)\ < Nc~ 2 . Let us show by contraposition that, with 
u) <£ B x>y (6), 

{ E : |G BL (x,y;£)| > 2a}n^(c) = 0. 

Assume otherwise, pick any point A* in the non-empty set in the LHS, and let 
J := {£" : \E' - X*\ < b} C M> C @(c). Then for any E 6 J one has, by (73J), 

> |/(A*)| - \J\ sup \f'(E')\ >2a- Nc- 2 ■ b > a, 
E'eJ 

so J C <? x , y (a) and mes(^ Xjy (a)) > mes(J) = 2b > b, contrary to the choice of uj. 
Since the set 3%{c) is independent of y, the assertion follows from (|7.6p . □ 
Taking into account the FEMSA bound obtained in Section Vi k = Pk < 
L k i one can set, e.g., with a = I and k > 6Ad, for L G {L^, fc > 0}: 

a(Lk) = L- k ^\ b(L k ) = L~^\ c(X fe ) = L^*"*)^. (7.8) 
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Since k > 6Nd, one has c(L k ) < L k 30 l i ' . 

Theorem 7.2. If the property (Wl) holds true, then for any pair of separable balls 
Bi t (x), Bi fc (y), k > 0, and a bounded interval Jcl, one has 

F{3EeI: min{F x (£),F y (£;)} > } < C £ -fe(i+«) fc +(™+i>* 

Proof. Introduce the event B b (x), b = b(L).), relative to the ball Bi t (x) and defined 
as in Theorem 1 7. 11 similarly, define the event Bb(y) relative to the ball Bi fc (y) and 
set B b = B b (x) U B b (y). Further, let <? x = {E e I : F x > 2a(L fe )}, <? y = {E e I : 
F y > 2a(Lfc)}, as in Theorem 1 7. II For any u £ B b , 

dist (4,S(H Bli (x))) < 2c(L fe ), dist (<f y ,£(H Blfe ( y ))) < 2c(L fe )- 

Therefore, applying Theorem 15.31 we obtain: 
P{(f x n^ y ^ 0} =P{dist(# x ,^ y ) =0} 

< P { B b } + P { dist (s(H Bifc ( X) ), S(H Bifc (y) )) < 4c(L fe ) } 

<L-* (1+e) V Const 4 2Ar+1)d (c(L fc )) 5 
<Consti fc -- (1+e)fc+(2JV+1)d . □ 

(7.9) 

7.2. Spectral localization. The assertion of Theorem 17.21 has a structure similar 
to that of the MSA bound from the work by von Dreifus and Klein [19] • More 
precisely, it guarantees a decay rate of Green functions slower than exponential, 
but faster than any power- law. The main difference is that one can rule out, 
with high probability, only the pairs of singular balls which are separable (and 
not just disjoint). It is not difficult to adapt the well-known argument from |19] 
and prove that with probability one, all polynomially bounded solutions to the 
equation H(lu)iP = Etp are in fact square-summable. The latter property requires a 
Shnol-Simon type result on spectrally a.e. polynomial boundedness of generalized 
eigenfunctions. It will follow independently by RAGE (Ruelle-Amrein-Georgescu- 
Enss) theorems (see a detailed discussion along with a bibliography, e.g., in [T7]) 
from the dynamical localization proven in Section [9j 

8. From fixed to variable energy: Second approach 

8.1. The spectral reduction. The next result is a notational adaptation of The- 
orem 7.1 from [15] . formulated there for single-particle Hamiltonians. 

Theorem 8.1. let be given a ball B/,(x), L > 1, and numbers a{L), b(L), c(L), 
Vl > obeying (|7.3[) and such that, for some interval I , all E G /, 

P{F X (£) >a}<V L . (8.1) 

Set K = |Bi(x)|. Then the following properties hold true: 

(A) For any b > Vl there exists an event B b such that ¥{B b } < b~ x Vh and for 
any uj £ B v the set of energies 

Ma) = ^x(a; u) = {F X {E) > a} n I 

is covered by K < 3K 2 intervals Ji = [E~ , E^\, of total length |«7j| < b. 
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(B) The endpoints E^ art determined by the functions E i-> (l x \(H. Sh r u \—E) 1 l y ) 
in such a way that, for the one-parameter family A(t) := H Bi ( u ) +t 1, the end- 
points E i (t) for the operators A(i) (replacing Hb 1 ( u )/' have the form 

E±{t)=Ef{0)+t, * G R. 

Proof. (A) Fix a point y € <9~Bi(x) and consider the rational function 



^ y Aj — E ^— ' Aj — E 

i—i %—\ 

Its derivative has the form 

and has < 2K — 2 zeros and < K poles, so f y has < 3K — 1 intervals of monotonic- 
ity ii ; y, and the total number of monotonicity intervals of all functions {/ y ,y G 
a-B L (x)} is bounded by K < \d~B L (x)\(3K - 1) < \B L (x)\(3K - 1) < 3Jf 2 , so 

U ye9 - BL{X) {E : / y (£) > a} = uf =1 J h J t = [£r,£+] C I, 

where, obviously, J2i 1^1 ^ mes{£' : F x (iJ) > a}. 

(B) Consider a one- parameter operator family A(t) = H Bi ( u )(cj) + 1 1 . All these 
operators share common eigenvectors; the latter determine the coefficients Hi, so 
one can choose eigenfunctions ipiit) constant in t and obtain Ki(t) = Ki(0). The 
eigenvalues of operators A(t) have the form Xi(t) = Xi(0) + t. We conclude that 
the Green functions, with fixed x and y, have the form f x . y (E; t) = f x . y (E — t; 0), 
so that the intervals Ji(t) have indeed the form J;(i) = [E~ + t,E^ +t]. □ 

Theorem 8.2. Assume the property (W2). Let be given numbers a > and Vl £ 
[0, 1] smc/i that (|8.ip ZioWs irue /or aZi x G .Z-^, i/iera for any pair of weakly separable 
(e.g., distant) balls Bi(x), B^(y) and for any b > Vl, one has 

¥ {3 E E I : min{F x (^),F y (£;)} > a} < 2b~ 1 V L + h L (Ab), (8.2) 

where 

h L {s) :=K 2 CL A s B +C'L A 's B \ K = max{|B L (x)|, |B z (y)|}. 

Proof. Fix b > (clearly, b < Vl gives rise to a trivial bound P{ } < 1). Define 
the event £>f,(x) relative to the ball B^(x) as in Theorem 18. II Similarly, define the 
event Bb{y) relative to the ball B^(y), and let Bb = Bb(x) U£>b(y). Further, let 
<?x,y be the event figuring in the LHS of (|8.2p and note that 

F { <S x , y } < P { B b } + P { s x , y n B c b } < 2b~ 1 V L + P { s x , y n B c b } . 

It remains to assess P{<S Xjy fl 23£ }. By Definition 15.21 there is a ball A C Z of 
diameter R < 2NL and index subsets J\,Ji C [1,N] such that n\ := J7i > 
\Ji | =: n 2 and 

n^Bf'lxjun^Bf^y) CA, 

n Ji 4 JV) (x)nA = 0, (8.3) 

n^Bf ) (y)nA = 0. 
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Consider the random variables 

^a(w):=\A\- 1 ^2v{z]u), rj,(w):=V(z;w)-&(w), z e A. 
zeA 

Denote the tr-algebra generated by the random variables {?y 2 , z € A; V(u, •), u £ A} 
by 3^, a- Consider the conditional probability distribution function 

F u (t\$ vA ) = F{t A <t\$ vA } 

and its continuity modulus 

^a( s I&7,a) = sup ess sup (-F ?A (i + s I^.a) - F^ A (t | S^.a)) ■ 
teK 

Owing to the assumption (W2), for some C, C", A, A' , B, S' g (0, +00) we have 

Vae[0,l] ¥{is u (s\$ ri , A )>CL A s B } <C"L A 's B '. (8.4) 

Using the representation V(z;u>) = ^\(u>) 1 +rj z (to) in the set A, introduce the 
respective operator decomposition 

H Bi (x)H=nia(w)l+A x (w) 1 (8.5) 

where the operator A x (w) is non-random, conditional on Sij a- F° r an y w ^ ^b( x )j 
the energies £ where F X (.E) > a are covered by intervals Ji yX , with £\ |J»,x| < &• 
Combining (|8.5|) and the assertion (B) of Theorem 18. 11 we can write 

J»,x(w) = [J5r x (a;) + ni£ A (cj), ££ x (w) + ni£ A x(cj)] 

where E^ x (uj) are Sv^-measurable. For any cj ^ £>t,(y), the energies _E where 
Fy(-E) > a are covered by intervals Jj ; y, also obeying | J^ y | < b. As above, 

Ji, y (w) - [E^(u) + n 2 U(u), E+ y {u) + rc 2 £ A M] 

where Ef y (uj) are also S^.A-measurable. 

Let = Sft(x) U£>;,(y) and set e iiX = |Jj, x |, e Jiy = Jj, y |. For u e i3g , one has 
< e ilX , £i !y < 26, so 

: J 4 , x n Jj, y # 0} n Bl C { |Sr x ( w ) - ET y (u) I < e ilX + e jt y } D Bg 
C ||(m -na)^A(w) -^j(w)| < 26 + 2&J 



where Eij(u)) is 3,,^-mcasurable. Set n = ni — «2(> 1). By 

p{|t»fA(w)-f ( j(w)| <4&} =E[p{|n£ A (w)-i3u(w)| <4&|&,, A } 

< P { i/ a (46 I ^, A ) > CX A (46) B } + CX A (46) B 

Observe that the first probability in the RHS does not depend upon i and j, so it 
suffices to count it only once. Taking the sum over we obtain the required 

bound: 

p { s x , y n B c } < 2b~ 1 v L + P { u ■ ( Ji, x n j jtY ) ^ } 

<2b- 1 V L + h L {Ab). □ 



In particular, taking into account Theorem 16.61 we can set, for L = 

a = a(L k ) = e -y{m,L k )L^ = L --a+e)\ b = h{Lk) = L -*W 
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These settings give rise to the following corollary of Theorem 18.21 

Theorem 8.3. Assume the property (W2). If there is an integer Lq > 1 and 
numbers m > 1, a € (1, 2) such that 

min{^ ,Qo} < Cf N Lo K , k > 4§^, 

then for some C, a,c > and all k > 0, any interval I C M with \I\ < 1 and any 
pair of distant balls Bi t (x), Bi fc (y), the following bound holds true: 

¥{EeI: B Lfc (x) andB Lk (y) are (E,m)-S } < Ce"" ll,1+Cit . (8.6) 

8.2. On the validity of the assumption (|8 .4(1 . First of all, recall that, by an 
elementary result on Gaussian distributions, if V : Z x fl — > K is an IID Gaussian 
field, say, with zero mean and unit variance, the sample average £ x of the sample 
{V{z;uj),z e Bl(x)} is independent of the sigma-algebra generated by the "fluc- 
tuations" rj z (uj); moreover, it has Gaussian distribution 7V(0, |Bl(x)|) and admits 
a probability density with Hj^Hoo 

< -J= IBl(x)! 1 / 2 . In this particular case, Eqn 
(|8.4p can be replaced by a stronger, deterministic bound: the conditional continuity 
modulus !/£ x (s | § x ) is actually independent of the condition and is a.s. bounded by 
IbeJIoo ■ s. 

Such a situation is rather exceptional, as shows the example of two IID random 
variables Vi(uj),V2{oj) with uniform distribution Unif([0,l]). Indeed, in this case 
£ := (V\ +V2V2, r\ — {V\ — V2)/2 and the distribution of £ conditional on 77 is uniform 
on the interval I v of length 0(1 — [77D, hence, with constant density 0{ 1— \r)\\ ), for 
?7| < 1; for r\ = ±1, this distribution is concentrated on a single point. However, this 
example shows also how such a difficulty can be bypassed: excessively "singular" 
conditional distrubutions of the sample mean £ occur only for a set of conditions 
having a small probability. Using this simple idea, Gaume |27j , in the framework of 
his PhD project, established the property (|8.4[) for IID random fields with piecewise 
constant marginal probability density. By standard approximation arguments, it can 
be easily extended to piecewise Lipshitz (or Holder) continuous densities, which is 
sufficient for most physically relevant applications. We believe that some variant of 
the property (I8.4[) . perhaps weaker but still sufficient for the purposes of the MSA, 
holds true in a larger class of IID random fields. 

8.3. Exponential spectral localization. The assertion of Theorem 18.31 has the 

same form as in the conventional MSA bound going back to the work by von Dreifus 
and Klein [TO] (actually, even slightly stronger); therefore, the same argument as 
in [TO] (having its roots in 26 ) applies and proves that with probability one, all 
polynomially bounded solutions to the equation H(uj)ip = Ei\) are in fact decaying 
exponentially fast at infinity, thus the operator H(uS) has a.s. pure point spectrum. 
The latter property requires a Shnol-Simon type result on spectrally a.e. polynomial 
boundedness of generalized eigenfunctions; it will also follow by RAGE theorems 
from the dynamical localization proven in Section [5] 

9. From MSA to the strong dynamical localization 

This section is a straightforward adaptation of Section 8 from 15 . Formulated 
in the most general form, the results of this section apply indifferently to single- 
and multi-particle Hamiltonians, and the particularity of the latter resides in the 
geometry of pairs of finite balls involved. 
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Recall that the first rigorous derivations of the dynamical localization from MSA- 
type probabilistic bounds on the Green functions have been obtained by Gcrminct 
De Bievre [25] and Damanik-Stollmann [5D]. Germinet and Klein proposed 
a shorter proof, which can be further simplified in the context of finite-volume 
operators (which merely should have compact resolvent), and essentially reduced to 
an elementary application of the Bessel inequality. 

In applications to iV-particle systems on a graph obeying the growth condition 
(|2.1[) . the parameter D below has to be replaced by Nd. 

9.1. EF correlators in finite balls. Given an interval J CR, denote by 
the set of all Borel functions : R — s> C with supp <f> C I and ||^||oo < 1- 

Theorem 9.1. Fix an integer L G N* and assume that the following bound holds 
for any pair of disjoint balls Bl{x), B^ly) o-nd some quantity s{L) > 0: 

P{3E G I : B L {x) and B L {y) are (E, m)-S } < ?(X). 

Then for any x,y S Z with d(x, y) > 2L + 1 , any finite connected subgraph (of Z) 
Q D Bl(x) U Bi(y) and any Borel function <p £ 

E[\{l x \<p(H g (uj))\l y )\] <4e- mL + q (L). (9.1) 

Proof. Fix points x,y G Z with d(x, y) > 2L + 1 and a graph Q D Bl{x) U Bi(y). 
The operator Hg(uj) has a finite orthonormal eigenbasis {ipi} with respective eigen- 
values {Xi}- Let S = dB^ix) U dB^iy) (recall: this is a set of pairs (u, w')); 
note that |S| < 2CpL D , by (I2.1[) . Suppose that for some ui, for each i there 
is Zi G {x, y} such that B^(zi) is (A^, m)-NS; let {vi} — {x, y} \ {z{\. Denote 
Hx.y(4>) = | {la; \4>{Hg{u)))\ l y }|, with fJL Xt y{4>) < 1. Then by the GRI for the eigen- 
functions, 

Hx,y(4>) < ||<£||oo ^2 \^i(x)lpi(y)\ < ^2 \^i{Zi)lpi(Vi)\ 

<e-™ L Y,(c 2 D L D )-^ \Mv)\(M*)\ + \My)\) 

\i£l (,u,u')eS 
D \i£l 

(using Bessel's inequality and |S| < 2C\,LP) 

< e- mL 2 max (2|| l u f + || l x || 2 + || 1, f) = ^ mL . 

Denote «S L = {3E el : B L (x) and B L (y) are (£,m)-S}, withP{S L } < by 
assumption. Now we conclude: 
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9.2. Dynamical localization on the entire graph. For the reader's conve- 
nience, we repeat below a simple argument, described in [15] and employed earlier 
by Aizenman et al. [2H]. The quantities fJ>x,y {$) = (l^ \<t>(H) \ l y ) defined, for 
example, for bounded continuous or Borel functions 4>, generate signed (i.e., not 
necessarily positive) spectral measures associated with a self-adjoint operator H: 

df, x H J(E)<P(E) := (l x \4>(H)\ l y ). 

In particular, we can consider, with x,y,u G Z hxed, measures related to 
operators Hb l ( u )j for all k > 0, as well as their counterparts fi x ,y for the oper- 
ator H on the entire graph Z. A sufficient condition for the vague convergence 
Hxy^r ^x,y as k — > oo is the strong resolvent convergence H-q l (u) ~~ ► H. Such 
convergence is well-known to occur for a very large class of operators, including 
(unbounded) Schrodinger operators in Euclidean spaces and their analogs on the 
so-called quantum graphs. Indeed, for (not necessarily bounded) operators H n with 
a common core T> to converge to an operator H with the same core, it suffices that 
H^ip Hip strongly for any element rp G T> (cf. [32]). For finite- volume operators, 
one can usually find an appropriate core T> formed by compactly supported func- 
tions ip; for finite-difference Hamiltonians on graphs (even unbounded, e.g., for DSO 
with unbounded potentials) one can choose as T> the subset of all functions with 
finite supports. On such functions, Hb l (u)V' ~~ * Hip as k — > oo (by stabilization), 
therefore, the spectral measures converge vaguely: — > \i x ,y By Fatou lemma, 
for any bounded Borel set 4cl, one has 

\H X , V {A)\ <liminf \nt, y {A)\ 

(here |/x(j4)| := sup{/i(0), ||0|| < 1, supp <p C A}). Taking the expectation and 
using the uniform upper bounds on EF correlators in finite balls, we conclude that 



E 



sup \(l x \<t>(H(u))\l y }\ 



< Ce -aln 1+a d(x,y) (g^ 



(using the inequality L^ 1 ^ < Ce~ alnl+C L « , for some C, a, c > 0). In particular, 
with functions <p t '■ A i-> e _l * A , we obtain the strong dynamical localization property 
for the ensemble of random Hamiltonians H(u>). 

Taking into account Theorem 18.61 we come to the following sufficient conditions 
of strong ./V-particle dynamical localization in an interval ICR: 

Theorem 9.2. Assume (Wl) and suppose that, for some N > 2 and all k > 0, 
S(I, N, k) holds true. Then there are constants a, c G (0, -f-oo) such that for all 
x G Z N , some C(x) G (0, +oo) and all y G Z N 



E 



sup \(l x \e- imlN) ^Pi(HW(u))\l y )\ < C7(x)c- alnl+Cp(x ^. (9.3) 



Here and below, Pj (lj)) stands for the spectral projection of operator 
HW(u;) on the interval /. Recall that, formally, we assumed in Sections 0- 
12] the interval / to be finite. In the case of a bounded random potential, the 
spectrum of the operator is covered by a finite, non-random interval /, so that 
P r (H( JV )(u;)) = HW(a;). 
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Theorem 9.3. Assume (W2) and suppose that, for some N > 2 and all k > 0, 
S(I, N, k) holds true. Then there are constants C,a,c G (0, +oo) such that for all 



E 



supKl.le-^^PKHWM)!^}) 
10. Appendix 



<Ce~ aln Ps(^y). (9.4) 



10.1. Proof of Lemma 16.41 The ball B := Bi t (u) is assumed PI, so there is a 
decomposition of the configuration u into two non-interacting sub-configurations, 
u = (u',u"), so that U(u) = U(u') + U(u"), and the operator reads as 

follows: 

H b<>) = H <V) ® 1<""> + 1<»'> ® H <)(x „, (10.1) 
thus its eigenvalues ax 6 the sums E a h — A a + fib, where {A a } = E(H B (u ,)) is the 

spectrum of Hb l ( u ') an d, respectively, {fib} = S(Hb l ( U "))- Eigenvectors of ^ 
can be chosen in the form \P a & = (p a ®ip b where {cj) a } are eigenvectors of H „(»'), M 

and {V'b} are eigenvectors of H („«). For each pair (A a ,/i(,), the non-resonance 

B Lk (u ) 

assumption |.E — (A a + fib)\ > e Lk reads as \{E — X a ) — fib)\ > e L fc and also as 
\(E - fi b ) - X a )\ > e L k. Set G(E) = (H b (jv), - - E)~ l ; then we can write 

G(u , yiE) . £E *- (,;i t^-?' w m 

= ^P:(u',y')GB ItK) (u" iy ";i?-A a ) (10.3) 

= J2PU^",y")G BLk(un (u',y';E- fi b ), (10.4) 
Mb 

By assumption, for all fib G <t(H^ f u »0, the projection ball Bl^ (u') is m)-NS, 
and for all A a S afHg^ ( U '))' ^ e P r °j ec tion ball Bi fe (u") is (A a ,m)-NS. 

For any y G Q^Bj^u), either p(u',y') = Lk, in which case we infer from (|10.4j) . 
combined with (fib, m)-NS property of the ball Bi fc (u'), that 

|G(u,y;B)| < |B Lfc (u")|e- 7(m ' ifc ^- 1)Lfe+2 ^ (10.5) 
or p(u",y") = Lk, and then we have by (|10.3I) 

|G(u,y;£)| < |B ifc (u')| e -^ m ' Lk ' N -^ Lk+2L ^ . 
In either case, the LHS is bounded by 

exp (-m(l + L^f-^-V+iLu + 21% + Const lni fe ) < \ e -^ m ' L *> N \ 
for Lq large enough, since m > 1 and 

L^ 7 " = l£ /8 » L^ 2 + Const In L k = L{ + Const In . 

□ 
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